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Self-gravitating accretion disk in Sgr A* few million years ago: 
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Abstract. Sgr A* is extra-ordinarily dim in all wavelengths requiring a very low accretion rate at the present time. 
However, at a radial distance of a fraction of a parsec from Sgr A*, two rings populated by young massive stars 
suggest a recent burst of star formation in a rather hostile environment. Here we explore two ways of creating 
such young stellar rings with a gaseous accretion disk: by self-gravity in a massive disk, and by capturing "old" 
low mass stars and growing them via gas accretion in a disk. The minimum disk mass is above 10* Mq for the first 
mechanism and is few tens times larger for the second one. The observed relatively small velocity dispersion of 
the stars rules out disks more massive than around 10^ Mq: heavier stellar or gas disks would warp each other by 
orbital precession in an axisymmetric potential too strongly. The capture of "old" stars by a disk is thus unlikely 
as the origin of the young stellar disks. The absence of a massive nuclear gas disk in Sgr A* now implies that 
the disk was either accreted by the SMBH, which would then imply almost a quasar-like luminosity for Sgr A*, 
or was consumed in the star formation episode. The latter possibility appears to be more likely on theoretical 
grounds. We also consider whether accretion disk plane changes, expected to occur due to fluctuations in the 
angular momentum of gas infalling into the central parsec of a galaxy, would dislodge the embedded stars from 
the disk midplane. We find that the stars leave the disk midplane only if the disk orientation changes on time 
scales much shorter than the disk viscous time. 



1. Introduction 

The complex chain of events leading to the growth 
of supermassive black holes (SMBHs) in the galactic 
centers is not yet fully understood (e.g., iRee s 2002). 
Nevertheless, gas accretion is probably the dominant 
physical process delivering the food to the giant black 
holes (e.g., .Yu fc Tr cmainc 2002). On small radial dis- 
tances from the SMBH, th e standard thin accretion disk 
ijShakura fc SunvaevI IiOtI appears to be a sure way to 
provide the SMBH with gas at rates approaching the 
Eddington limit. However, at distances larger than ~ 10~^ 
parsec from the SMBH, standard thin disks run into sev- 
eral problems. 

First of all, the time necessary for the gas to inflow 
into the black hole - the disk viscous time scale - be- 
comes too long (e.g. as large as 10'* — 10^° years for 
larger radii). In parallel with this, the (standard accretion) 
disk mass becomes very large. When the latter exceeds 
about 1 % of the SMBH mass, local gravitational instabil- 
ities d e velop (e.g..|Paczvnskll978l:fKqlvkhalov fc SunvaevI 



i.g. 

1980' Shlosman 



Beeelmanl Il989l: iGoodmanI l20o' 



Collin fc Zahn.. 1999,1 . The structure of such disks is very 
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much open to discussion (see ^ due to uncertainties in 
theory and a dearth of relevant observations. 

Sgr A* is the closest SMBH (wi th a mass Mbh ~ 
3 X lO^M©; e.g.. lSchodel et al.ll2002^ . Although Sgr A* 
remained extremely dim during the entire history of X- 
ray observations, there are hints that it was much more 
active in the past. X-ray and 7-ray spectrum of the gi- 
ant molecular cloud Sgr B2 is most naturally explained 
as a time-delayed re flection of a source w i th a flat AGN- 
like spectrum (e.g.. ISunvaev et all Il993t iKovama et alJ 
Il996t iRevnivtsev et alJ l2004|) . The required luminosity 
is in the range of ~few xlO^^ erg/sec, too high by the 
Galactic standards. Sgr A*is then strongly suspected of 
being brighter in X-rays by some 6 orders of magnitude 
300-400 years ago. 

Deeper in the past, some few million years ago, few 
dozens of massive stars were formed and are currentl y at 
a distance of order 0.1 — 0.3 pc from Sgr A * (|Krabbc e t alJ 
ll995HGenzel et al.l2003llGhez et al.l2003l) . This is surpris- 
ing given that the tidal force of the SMBH would easily 
shear even gas clouds with densities orders of magnitude 
higher than the highest density cores of GMCs observed 
in the Galaxy. 

In situ star formation scenarios for the Sgr A* youn; 
massive stars have been numerically studied by ISander; 
l|l998l) with a sticky-particle code and also qualitatively 
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described bv lMorris. Ghez. fc Beckim ('l999"). In particu- 
lar, one of the simulations done by Sanders ( 1998 ) assumed 
that a cold cloud of gas with radial dimensions of 0.4 pc 
and with a small angular momentum infalls into Sgr A* 
gravitational potential starting from distance of 2.4 par- 
sec. The cloud gets tidally sheared into a thin band of 
gas which then forms a precessing eccentric ring. Frequent 
shocks are assumed to lead to strong compression of the 
gas and star formation. 

The existence of such low angular momentum clouds 
seems to be in question. In addition, the initial condi- 
tions of the simulations are rat her extreme: t o be stable 
against tidal shear (e.g., eq. 1 of lSanderlll998|) . the cloud 
mass should be Md ^ 7 x 10^ Mq, a very larg e mass for a 
cloud of 0.4 pc in size. The recen t obser vations ( Lisztl2003l) 
seem to contradict the 'Sanders' model for the ion- 

ized gas streamers. iGenzel et al. (2003) discount a current 
star formation in the mini-spiral, which is believed to be 
an ionized streamer. iGenzel et al.l l|200 t^ also note more 
generally that "if massive stars are forming frequently in 
dense gas streamers when outside the central parsec and 
then rapidly move through the central region, one would 
expect ~ 100 times as many massive stars outside the cen- 
tral region as in the central parsec" , which is not the case 
observationally. 

Al ternatively to the in situ star formation, iGerhardI 
ll200l|) proposed that the young massive stars could 
have been formed outside the central parsec in a mas- 
sive star cluster. Then, due to dynamical friction with 
the older population of background stars, the cluster 
would have been dragged into the central parsec and 
then dissolved there by the SMBH tidal shear. However , 
this appears to be only possible ijKim fc Morrid l2003t 
iMcMillan fc Portegies Zwar3l2003() if the cluster is very 
massive (M ^ 10^ M©), or if it is formed very near the 
central parsec already. In both cases a very dense core 
for the star cluster is required and appears to be unre- 
alistic. An intermediate mass black hole in the center of 
the cluster does allow the star cluster to survive longer 
against tidal disruption and hence transport the young 
st ars in the central parsec more efficiently (as suggested 
bv lHansen fc Mil osavlievicll2003l). However the numerical 
simulations of^l ^^^Figer. fc MorrisI l|2004l) show that the 
mass of the black hole has to be unusually large (~ 10%) 
compared with the cluster mass for this idea to work in 
practice; 

iLevin fc Beloborodovl l)2003|) and Nayakshin, Cuadra, 
fc Sunyaev (2004) suggested that the origin of the young 
stars is a massive self-gravitating accretion disk existing 
in Sgr A* in the past. Here we intend to investigate this 
idea quantitatively and to also look into some related the- 
oretical questions. 

We first estimate the minimum mass of such accretion 
disks to be around 10^, for each of the two stellar rings. 
In addition, we rule out the possibility that a less mas- 
sive accretion disk could capture enough of low mass stars 
from the pre-existing "relaxed" Sgr A* cusp and then grow 
them by accretion into massive stars {^^. 



We then attempt to understand the spatial distribu- 
tion of the young stars. In particular, we find that the rate 
of A^-body scattering between the stars (21) of the same 
ring can explain the observed stellar velocity dispersions 
in the inner stellar ring if the time-averaged total stellar 
mass in the ring was 10** Mq or higher. We also find that 
stellar orbits in both rings should remain close to the cir- 
cular Keplerian orbits up to this day (if stars were indeed 
born in a disk). The outer ring is however observed to be 
geometrically thicker and with a higher velocity dispersion 
than the inner one. The velocity dispersion of the outer 
ring may result from the stellar disk warping in the gravi- 
tational potential of the inner ring. Such warping sets the 
upper limit on the disk mass of about 10^ M© (^. 

We also question in ^ whether it is possible for the 
disk to leave the newly born stars behind (due to their 
high inertia) when the disk plane rotates. Dislodging the 
newly born stars, or proto-stars, from the disk midplane 
would have significantly reduced the problems faced by 
accretion disks at large radii since these stars would then 
stop devouring the disk and instead heat it and speed up 
the accretion of gas onto the SMBH via star-disk collisions 
(Ostriker....l983) . However, we find that the disk maintains 
a firm grip on these stars unless the plane change occurs 
on a time scale much shorter than the disk viscous time 

It is found that young massive stars would not mi- 
grate much radially ([Q in the disk, meaning that they 
are probably located at the radius where they were orig- 
inally formed. Small scale proto-stellar disks around the 
embedded stars may be gravitationally unstable as well 
and may create further generations of stars. Hierarchical 
growth and merging of such objects may result in the cre- 
ation of "mini star clusters" with the central object col- 
lapsing to an intermediate mass black hole (SQ. This could 
potentiall y be relevant to the o bservations of such objects 
as IRS13 (iMaillard et al.ll2nn4l) . 

Since the combined mass of the stellar material in the 
obser ved stellar rings presently is ^10'^ Mq (jGenzel et alJ 
there is then an interesting question of whether 
most of the gaseous disk mass has been used to acti- 
vate the presently dormant Sgr A* or it was reprocessed 
through star formation and expelled to larger radii via 
winds and supernova explosions. We believe the latter out- 
come is more likely since the accretion of gas onto embed- 
ded stars is very efficient. We briefly discuss observations 
that could distinguish between the quasar and the nuclear 
starburst possibilities (fQ. 

2. The minimum mass of a self-gravitating disk in 

Sgr A* is 10^ Mq 

The s tandard accretion disk solution l)Shakura fc SunvaevI 
Il973l) neglects self-gravity of the disk. Clearly, this solu- 
tion becomes invalid when the disk becomes strongly self- 
gravitating, but here we only want to estimate the mini- 
mum disk mass at which the self-gravity becomes impor- 
tant. For numerical values of the standard disk param- 
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eters, we follow ISvensson &: Zdziarskil l)l994|) with their 
parameters C = 1 ^nd / = (i.e. no X-ray emitting disk 
corona is assumed). The dimensionless accretion rate m is 
defined as m = M / M-^dd, where M-^dd = ^Edd/^c^ is the 
accretion rate corresponding to the Eddington luminosity 
and e ~ 0.06 is the radiative efficiency of the standard 
accretion flow. ^ 

For large radii [r 1) the gas dominated equations 
are appropriate: 



BH mass = 3x10' 



-1/10^1/20-1/5 



R 



n-2 ^-4/5 A,fy5^-3/5^3/5 

O 



T = 6.3 X 10^ K (aMs)-!/^™^/^ 



R 



-1 -9/10 



lO^Rs 



(1) 
(2) 

(3) 



Where H is the disk vertical height scale, R is the dis- 
tance from the SMBH, T is the midplane gas tempera- 
ture, a is the dimensionless viscosity parameter, Ms = 
Mbh/10*Mo, r = R/Rs, Rs = 2GMbii/c^ is the 
Schwarzschild radius of the SMBH and S is the sur- 
face density of the accretion disk. These equations as- 
sume Thomson electron scattering opacity for simplicity. 
Figure ^ shows some of the disk parameters obtained for 
-^BH = 3 X 10^ Mq, a — 1, and two values of the dimen- 
sionless accretion rate, to — 0.03 (thick lines) and to = 1 
(thin lines). The mass of the disk as a function of radius 
is approximated as Md ~ ttR^Y,. 

It is well known that the standard accretion disk be- 
comes self-gravitating at large radii, when the iToomrel 
l)l964|) parameter Q becomes less than unity, 



Q 



c^n H Mbh 



rGS R Md 



< 1 



(4) 



(cs is the sound speed inside the disk and fl its angular ve- 
locity) . The radius where Q = 1 yields the minimum mass 
of the disk needed for the latter to become self-gravitating. 
As can be seen from Fig.^ the disk should weigh at least 
10"* Mq in order to become sclf-gravitating. Note that this 
minimum disk mass estimate is quite robust because H/ R 
depends on a, radius and the accretion rate only weakly. 

This estimate is also conservative. The basic Shakura- 
Sunyaev model used here does not include irradiation by 
the central source, which may increase the disk midplane 
temperature somewhat, leading to a slightly larger H/R. 
In addition, trapping radiation by opacity effects would 
reduce the efficiency of cooling, adding to the stability of 
the disk against self-gravity. 

3. Capturing low mass stars and growing them by 
accretion: too slow. 

lArtvmowicz. Lin. & Wample 3 lll993h noted that stars in 
the nuclear star cluster on orbits relatively close to the 
local circular rotation of the accretion disk in quasars will 



Note that our definition of m is different from the one 
ISvensson fc Zdziarski il99'j) use: rhsz ~ 17.5m . 
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Fig. 1. The disk mass in units of 10^ M0, the Toomre 
instability parameter Q, and the ratio of the disk height 
scale H to radius i? as a function of radius for the standard 
accretion disk model. The thick curves are plotted for to = 
0.03, whereas thin ones are for m = 1. In both cases the 
disk is unstable to self-gravity at i? ^ 10^i?s, where its 
mass is Md > IO^Mq. 

be captured by the disk. The stars can then rapidly grow 
by gas accretion, and then enrich the accretion disk with 
heavy elements through stellar evolution. For the problem 
of the observed young massive stars in the Galactic Center 
(GC), the Artymowicz et al. trapping mechanism may be 
an alternative route to form the stars. The accretion disk 
does not have to be self-gravitating for the mechanism to 
work, provided there is enough stars and the time scale for 
the star trapping is right. One may thus hope to reduce 
the required disk mass. 

To within a factor order unity, equation 15 of 
lArtvmowicz et al yields the number of stars cap- 

tured by the disk within time At as 



/■4 

dN,{R) - ^N.,{R) 



(5) 



where A^* is the total number of stars in the star cluster 
within radius i?, and the variable ^ is defined by 

MMd At 



C = 32Cd- 



P 



(6) 



where P is the orbital period and M* = to* Mq is the 
typical mass of the stars in the cluster. For an estimate, 
we take Cd ~ 3 and q_2 = lOOMd/MsH ^ 1- Note also 
that we want to start with abundant stellar seeds, so we 
assume to* ^ 1. At the typical radial position of the young 
massive stars in Sgr A*, i? = O.lpc ^ 3 x 10^i?s, the 
circular Keplerian rotation period is 



P ~ 10^ year 



R 



-1 3/2 



3 X 105i?s 



(7) 
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Thus, 



At 



(8) 



Now, from results of lGenzel et alJ l)2003(l we estimate that 



N,(R) ~ 10^ ml 



R 



3 X 10^ Rs 



3/2 



(9) 



(we assumed the cusp power- la w index p = 1.5 for sim - 
plicity; see also equation 15 in iNavakshin et alJ l2004l) . 
Therefore, the number of stars captured by the disk is 



7 9_2 



At 



10^ years 



(10) 



Note that the radius R and the average stellar mass 
scaled out of this relation. 

The number of captured stars is somewhat low if we 
take into account the fact that the ~ 12 "Helium" stars 
found in each of the rings are only the brightest end 
of the stellar distribution, and there are pr obably many 
more (less massive) stars in the rings (§3.7 in lGenzel et alJ 
l)2003|) '). Therefore one would require q-2 ^3 — 10, that 
is disk mass ~ (1 - 3) x 10^ Mq. With this rather 
high required disk mass, the disk would have to be self- 
gravitating, and one expects a large number of stars being 
born inside the disk. The disk capture mechanism thus 
fails to reduce the minimum disk mass. Nevertheless, it 
should not be forgotten completely because of its ability 
to bring in some late type stars into the plane of the disk, 
e.g. the plane of the young massive stars. This type of stars 
would not be born inside the accretion disk in a time span 
of just a few million years. 

4. Velocity dispersion of an isolated stellar disk 

Stars embedded in accretio n disks are often considered 
in a "test star" regime fe.g.. lSver et alJll99l|) . when each 
star co-rotates with the accretion disk. The star's velocity 
is then nearly equal to the local Keplerian circular veloc- 
ity. When number of embedded stars, N^, ^ 1, two-body 
interactions between stars will increase their local veloc- 
ity dispersion, a, potentially leading to some interesting 
consequences. 

Since the disk velocity field has the radial Keplerian 
shear, it is the radial velocity dispersion of stars that 
will grow the fastest. However, when the anisotropy 
(Tr/<^z becomes larger than 3, the buckling instabil- 
ity will develop and the s tellar velocity dispersion 
will become mo re isotropic fe ulsrud. Mark, fc Carusol 
I1971I: IShlosmaiTfc Begelman 1989). Thus we assume an 
isotropic velocity dispersion here for simplicity. The veloc- 
ity dispersion of stars grows due to A^-body interactions 
at the rate 



da AnG^AUp^ 



In A, 



(11) 



dt (t2 

where A* ~ H^a"^ /GM^ is the Coulomb logarithm for 
stellar collisions; is the stellar disk height scale, which 



in general may be different from the gas disk height scale 
H. 

The growth of velocity dispersions opposed by the dy- 
namical friction force acting between the stars and the 
gas. Consider a star moving inside the disk with a relative 
vel ocity Vre} with respe ct to the local Keplerian velocity, 
wk- lArtvmowicd l)l994l) shows that the angular momen- 
tum and energy flow between the disk and the star (a 
small disk perturber in the case of proto-planetary disks) , 
calculated explicitly, coincides within a factor of few with 
the hydrodynamical Bondi-Hoyle drag acting on the star 
during its passage through the disk. The acceleration ex- 
perienced by the star is thus 



^)3/2 



l/4„ 



(12) 



is of the or- 



where g = min(l, Wc/wroi) and Vc = C, 
der of the escape velocity from the star, z;esc,*- In pertur- 
bative analytical approaches, such as dynamical friction, 
Cd 1 is the Coulomb logarithm. In A, where A is the ra- 
tio of the disk height scale H to Bondi (or accretion) radius 
('e.g.. lOstrikeJl983|) . However, in many circumstances the 
Bondi-Hoyle formula for accretion rate onto the star pro- 
duces supper-Eddington values. The drag force (e.g. Cd) 
should then be reduced to account for the radiation pres- 
sure force. One finds that in the disk geometry the largest 
contribution to the star-gas friction comes from distances 
^ H from the star. The exact value of Cd depends on disk 
opacity and the 3-D velocity of the star, but estimates 
suggest that Cd is not much smaller than unity in this 
case. 

Note that when the relative velocity is high, the drag 
force is just the hydrodynamical drag, Od oc irR^pv^^y, 
where R <t is the stellar radius. For Cs < Wrci < fcsc,*, the 
classical IChandrasekhaJ l)l943(l 's dynamical friction for- 
mula is recovered, with Od oc v~2- Finally, if the relative 
velocity is smaller than Cs, we have ~ MBWrcb which is 
about equal to the momentum flux accreted by the star 
(Mb is the Bondi accretion rate). 

While Drci is not too large, i.e., 3=1, the evolution of 
the stellar velocity dispersion is approximately given by 



dt 



: In A* 



pCdCr 



(c2 -1-^2)3/2 



(13) 



As long as p*lnA* < pCd, the star-gas drag will be 
able to keep the stars on local circular Keplerian orbits 
in the sense that ti <C Cs, the gas sound speed, thus 
the stars indeed behave as test particles. However, when 
/5H.lnA,t > pCd, the stellar velocity dispersion will evolve 
mainly under the influence of A^body collisions, and it will 
run away. 

It may appear that the last fact suggests a natural 
mechanism for stopping the very efficient (see Sj^J accre- 
tion of gas onto the embedded (proto-) stars. As the stel- 
lar velocity dispersion grows much larger than the gas 
sound speed, the stars will be no longer embedded in ac- 
cretion disks as they would spend most of their orbits 
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outside the main body of the accretion disk. In addition, 
even when the stars are crossing the disk, the relatively 
high value of frci means that the accretion rate onto stars 
will be strongly reduced. However, the effect is important 
only when > p (assuming In A* Cd), that is when 
the stellar density is already larger than the gas density. 
Therefore, before this effect may become important, about 
a half of the initial gas mass should already be consumed 
by the stars. The accretion onto the stars is curbed by 
the A^-body dispersion effects too late, when the disk is 
already half eaten by the stars. 

Now, coming back to the Sgr A* case, we can estimate 
the expected ~ Ra/v^. The relaxation time, defined 
as the time needed for the stellar disk to thicken to height 
iJ* can be found from equation 1111 



^dyn 



1 4 



4Afd,M*lnA, 



(14) 



where Ma* is 7iB?H^p^^ the mass of the stellar disk. 
Equation yields 



^rcl 
^dyn 



2500 



lO^Mp, 10 Mp, 



d* 



-ff* / R 



0.1 



InA: 



(15) 



With tdyii ~ 300 years, we have tj-ci ~ 10^ years. Hence 
the geometrical thickness of the rings, and the ratio of 
velocity dispersion to the local Keplerian velocity, (t/wk, 
are expected to be of order 0.1 for the two young stellar 
rings in the GC. The individual stellar velocities should 
thus be still close to the local circular Keplerian values if 
th e origin of the stars is in the gaseous disk. 

iLevin fc~B cloborodov ( 2003) estimate the geometrical 
thickness of the inner stellar disk in Sgr A* to be of or- 
der H^,/R = 0.1. This ratio is howe ver larger but is not 
quantified for the outer disk found bv lGenzel et alJ (|^03) . 
From their figure 15 we estimate that Hf,/R ^ 0.3 for the 
outer, counter-rotating, disk. 

One may try to invert equation 1151 to constrain the 
initial stellar mass of the dis ks in the GC by us ing the 
observed velocity dispersions l|Genzel et al.lEo03h in the 
rings. Unfortunately the limits are not very stringent due 
to the strong dependence of on H/ R. A disk mass as 
high as Md 3 x 10^ Mq could still be consistent with 
the observations for the inner stellar ring. Interestingly, for 
the outer stellar ring, the velocity dispersion is too high 
to be explained by the A^body effects unless the ring mass 
is unrealistically high. 

5. Destruction of stellar rings by orbital 
precession: the maximum disk mass. 

iGenzel et al.l ( 20031) find that most of the young innermost 
stars lie in one of two stellar rings. There is no noticeable 
difference in the estimated age of the two groups of stars. 
The rings are bound to interact gravitationally with one 
another, and this could lead to observable disk distortions. 

In particular, stellar orbits precess around the axis 
of symmetry in an axisymmetric potential (e.g., §3.2 in 



'Binn ev fc Tremaind Il987j) . We represented one of the 
disks by the Kuzmin potential 



^k{R,z) 



(16) 



v/i?2 + (a + |z|)2 

where a is the disk radius, R is the radius in the cylindrical 
coordinates and z is the perpendicular distance from the 
disk. We then numerically integrated stellar orbits, start- 
ing from nearly circular Keplerian orbits unperturbed by 
the disk presence. The orbits remain approximately circu- 
lar, and conserve the inclination angle i between the or- 
bital plane and the disk plane (because the z— component 
of the angular momentum is rigorously conserved in the 
axisymmetric potential). The stellar orbital plane pre- 
cesses with respect to the disk at a rate 

4) ^C^qP-^ cos i , (17) 

where Cp is a constant (for a given orbit and given geom- 
etry) of order unity. Angle (j) here is the azimuthal angle 
of the lines of the nodes for the orbit in cylindrical coordi- 
nates used to define the Kuzmin potential. The precession 
rate scaling f equation I17|l is natural since for small q the 
effect is linear in q as can be seen for orbits nearly co- 
planar with the disk (when i « 0°); for i = 90° there 
should be no plane precession due to symmetry. 

The value of Cp depends on the value of a with respect 
to the radius of the nearly circular stellar orbit; for a of 
order the radius, Cp ~ 1. S etting i = 74.° as a ppropriate 
for the two GC stellar rings (|Genzel et alJ2003|) . we obtain 
q t 



Acj)^ Cp 



(18) 



P 0.003 103P ■ 
The important point to note is that nearly circular 
orbits of stars at different radii from the SMBH will pre- 
cess by different amounts A(f>. Therefore such a precession 
leads to a warping of the stellar disk. After a time long 
enough to yield A(f) 1 somewhere in the disk, the ini- 
tially flat stellar disk will be disfigured and will not be 
recognizable as a disk at all by an observer. 

An approximate upper limit on the time-average mass 
of each of the two stellar rings in Sgr A* can be set. 
Clearly the exact value of this limit should be obtained nu- 
merically with A^body experiments and comparison with 
the quality (y^) of the fits to the two observe d planes 
llLevin fc Beloborodo-^l2003t ICenzel et al1l2003() . Such a 
study is underway. Due to an observational uncertainty in 
the radial dimensions of the rings' inner and outer radii, 
and theoretical uncertainty in the distribution of gaseous 
mass (i.e. S(i?)) in the accretion disk, it is possible to 
reduce Cp from its maximum value for some values of pa- 
rameters. Nevertheless, a rather robust value for the upper 
mass of the disks appears to be 

maxMd« 10^ Mq. (19) 

6. Rotating the accretion disk midplane: do stars 
remain embedded? 

The accretion disk midplane orientation can in principle 
change as a result of a new mass deposition coming with a 
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different orientation of the angular momentum vector, fn 
sucli a rotation, would the newly born stars remain embed- 
ded into the disk and follow its rotation or would they stay 
behind in the "old" accretion disk midplane due to their 
large inertia? The answer to this question is important for 
AGN disks in general as embedde d stars can seriously in- 
fluen ce the accretion process fe.g- lGoodman fc Tanll2004l: 
lNavakshin,2004) . 

By order of magnitude, one can estimate the time 
needed to turn the accretion disk plane around on a sig- 
nificant angle to be 



Id M 



i-visc • 

Mr, 



(20) 



where is the mass condensation rate onto the accretion 
disk. If the accretion and condensation processes are in an 
approximate steady state, — M, t^ot 
is 



tvisc- The latter 



^visc ^ 
VK 



(21) 



and can be fairly long. Thus in general the disk plane 
orientation changes rather slowly. 



6.2. Critical rotation time 

Suppose that the accretion disk midplane turns at a rate 
given by the time scale trot • Define a critical rotation time 
scale, ire, such that for disk plane changes occurring on 
time scales shorter than tic, the stars are dislodged from 
the gas disk. For tiot > trc, on the contrary, the stars re- 
main bound to the disk. Clearly, we get the critical time 
scale when Cpi ~ a^ot = I'^/trot, where a^ot is the "rota- 
tion acceleration" of the turning disk midplane. We obtain 
for the critical rotation time 



_ Mbh 
~ 2Md ' 



yn • 



(24) 



Figure |2] shows the critical rotation time scales (dotted 
curves) along with other important time scales for the 
standard accretion disk model with same parameters as 
used for Fig. and for a 10 Solar mass star. The thick line 
curves are for m — 0.03, whereas the thin curves are for 
TO = 1. The accretion and migration time scales will be 
discussed in 13 below. 

Note that trc is longer than tdyn but is much shorter 
than tvisc- This implies that if accretion disk plane changes 
occur on a viscous time scale, the stars would remain 
bound to the disk. Only very fast plane changes could 
dislodge the stars from the disk midplane. 



6.1. Forces l<eeping the stars embedded 

Two forces mediating interaction between a star and a 
gaseous disk are the gravity of the disk as a whole, and 
the friction acting on a star moving inside the disk at a cer- 
tain velocity with respect to the local circular Keplerian 
speed. If stars lag behind the rotating disk plane, the char- 
acteristic relative velocity at which the star and the gas 
would be separated is ~ WK/^iot ^^nd is very small com- 
pared to the sound speed in the gas (if trot ^ tvisc)- At 
small relative velocities the dynamical friction of a star 
"leaving" the gas disk is very small too (see eauation ll2l) . 
and a simple estimate shows that the dynamical friction 
force can be safely neglected in what follows below. 

Therefore the binding force to consider is the direct 
gravitational attraction between the disk and the star. 
Near the disk midplane, the infinite plane approximation 
can be used for the disk gravity. The gravitational attrac- 
tion force of the gaseous disk for a star that left the disk 
midplane (i.e., the star-disk midplane separation |2:| ^ H) 
is 



27rGS 



(22) 



where S is the local disk surface density. Comparing this 
acceleration with the centrifugal acceleration of the star 
moving in a circular Keplerian orbit around the central 
black hole, Oc — v^/R, 



dpi 



2Md 
M, 



BH 



(23) 



6.3. The case of Sgr A* 

We have just shown that it is fairly difficult to separate the 
stars and the accretion disks in slow disk plane rotations 
or deformations. For the Sgr A* case, this implies that 
either (i) there were two separate accretion disk creation 
events that created the two differently oriented rings; or 
(ii) the accretion disk itself was extremely warped so that 
its inner part was oriented almost at the right angle with 
respect to the outer disk part. 



7. Accretion onto embedded stars 

The Hill's radius Rh, 



Rh 



3AfBH 



1/3 



R, 



(25) 



defines the sphere around the star where the dynamics 
of gas is dominated by the star rather than the SMBH. 
The accretion of gas onto a star is believed to be simi- 
lar t o the growth o f te rres trial plane ts in a planetesimal 
disk l(LissaueillT987HBate et al.ll20o"'j) . For Rh > H, gas 
accretion onto the star is quasi two-dimensional. The ac- 
cretion rate is determined by the rate at which differential 
rotation brings the matter into the Hill's sphere, 



= Mh ^ AttRhH pvH ^ AttRhPCs , 



(26) 



where p — Y./2H is the mean disk density. We used 
the fact that the characteristic gas velocity (relative to 
the star) at the Hill's distance from the star, v^, is 
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Radius, Rg 

Fig. 2. Time scales for a 10 Solar mass star embedded in 
the standard accretion disk with parameters as in Fig. ^ 
As before, thick curves correspond to m = 0.03, whereas 
thin ones are for rh ~ 1. The solid lines show the viscous 
and the dynamical time scales for the disk, as labelled 
in the figure. The star is massive enough to open up a 
gap and hence migrates inward on the viscous timescale 



The dashed and dotted lines are the ac- 



cretion and the critical rotation time scales, respectively. 
For chosen parameters, the former one is independent of 
rn (see text in Sj3for detail). 



Vh = RH\dfl/d\nR\ ~ Cs{Rh / H) since the angular ve- 
locity for Keplerian rotation is = c^/iJ. Eauationl25lis 
valid as long as Rh > H since in the opposite case the 
gas thermal velocity becomes important and the accre- 
tion would proce ed at the Bondi accretion rate (Mb; e.g., 
ISver et ablFlim]) . Of course Af* cannot exceed AT^.Edd — 
10~^r, ~ lO~^M0/year mi^^, the Eddington accretion 



rate onto the star ^. We thus estimate 

= min \Mh, A>/b, A>/. Edd 



(27) 



One can then define the accretion time scale for a star 
embedded into a disk: 



(28) 



Figure 121 shows the accretion time scale (dashed line) for a 
10 Mq embedded star. Although we considered two values 
for the accretion rate onto the SMBH, m — 0.03 and m = 
1, as in Figure n tacc turns out to be the same for both of 
these because the accretion rate is close to the Eddington 
value. 



We assumed th at r<. = {Rtl Rc- 
that equation 10 in lNavakshinI i2004l 



) « (M./Mq)!/^ Note 
contains a typo. 



An important point to take from Figure [3 is that ac- 
cretion onto embedded stars is able to double the stellar 
mass in few thousand years. Therefore, growing stars as 
massive as 100 Mq in a million years in a disk with gas 
mass Md ^ 10* Mq appears to be no problem at all. [One 
potential uncertainty here is the reduction in the accre- 
tion rate onto the embedded stars once these stars are 
massive enough to clear out a radi al gap in the accretion 
disk. Results of iBate et all ()2003|) . Fi gure 9, show that 
this reduction can be very large. However, in the case of 
an AGN disk with many embedded stars, the dynamics of 
the star-gas interaction is surely going to be different from 
the case of a "test" star or planet. The accretion disk will 
then be divided onto many rings between stars on nearly 
circular radial orbits. If the orbits are close enough (num- 
ber of stars Ni, ^1), then the gas in a ring will experience 
alternating inward and outward pushes from the two stars 
closest to it and hence the radial gap can in fact be closed, 
enabling an unhindered accretion. The issue deserves a fu- 
ture study.] 

We also estimated the radial migration time scale, 
^migr, using thc prescription for the radial migr ation ve- 
locity based on the numerical calculations of Bate et alJ 
l|2003() . For the parameters chosen, the 10 Mq stars, and 
any stars more massive than that, open up a gap in the ac- 
cretion disk and their radial migration is identical with the 
viscous flow of matter in the gas disk. Thus tmigr = ^visc 
(two solid curves in Figure |2J|. The migration time scale 
is very long, indicating that stars will remain pretty much 
where they were born in accretion disks with parameters 
close to that of the standard disk for Sgr A*. In fact, a 
more realistic self-gravitating disk would not change this 
conclusion significantly since the migration time scale only 
gets longer when the midplane disk density decreases as a 
result of disk swelling due to gravitationally induced tur- 
bulence. 

8. Growth of "mini star clusters" and 

intermediate mass black holes in accretion 
disks 

iGoodman fc TaiJ (|20^ have recently suggested that it 
is possible to grow supermassive stars in AGN accretion 
disks. The maximum mass of a star in this case is the gas 
disk mass in a ring with width of order the Hill's radius 
of the star, Rh = (Af*/3A^BH)^/^- This is the "isolation" 
mass, Af, M^^'^Mr^''^ 



BH 



A/, « 550 M: 







Md 


3/2 


■3 X IO^Mq' 


1/2 


IQ'^Mq 




A/bh 





(29) 



The stability of super-massive stars is briefly summarized 
in §2 of (Goodman & Tan 2004). The supermassive star 
could collapse directly i nto a black hole i f the star is more 
massive than 300 Mq l)Frver et al.ll200l|) . 

However, the Hill's accretion rate estimate assumes 
that all of the disk mass delivered by the differential rota- 
tion into the Hill's radius about the star is accreted onto 
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the star. Even without the gap presence, this is not ob- 
vious because the gas still has to loose most of its an- 
gular in£nientum__befoi^^ the stellar surface 
(;e.g.. lMilosavlievic fc L'oei3l200 4') . Furthermore, quite fre- 
quently the accretion rate onto the star estimated in this 
way exceeds the Edding ton accretion rate onto the st ar (as 
is the case for Figure 2').'Milosavlievic fc Loebl l)2004') have 
shown that the fringes of the small scale disk around the 
embedded stars themselves become self-gravitating and 
may therefore also form stars or planets. It is thus possible 
to grow in situ star clusters. The maximum total mass of 
such a cluster should be close to the isolation mass. 

A qualitative confirmation of these ideas can be found 
in numerical simula tions of a related physical problem by 
iTanga et al.l l|20fl4fl . These authors simulate the growth 
and clustering of planetesimals in a proto-stellar disk. 
They find a hierarchical growth of clusters of particles 
and find that these "clusters" are intrinsically stable struc- 
tures. This is likely because of the abundant supply of gas 
into the Hill's sphere: there is always a plenty of gas to 
interact with the particles (gravitationally bound objects 
in AGN case) that continue to get more and more bound 
to the central object in the mini-cluster. 

This mechanism of intermediate mass black hole 
(IMBH) and bound to it star cluster may be rele- 
vant t o the observat i ons o f the IRS13 cluster near 
Sgr A* l|Maillard et al.ll2004l) . The "dark" mass in the 
IRS13 is estimated to be > IO^Mq. Equation shows 
that an initial mass of the disk of order several times 
10* Mq would have been sufficient to grow "in situ" an 
object massive enough to become the IRS13 cluster. 

9. Discussion 

We have considered here the formation of massive stars in 
a self-gravitating accretion disk for conditions appropriate 
for the central ~ 0.2 pc of our Galaxy. Formation of an 
accretion disk (instead of a narrow ring) would be a likely 
outcome of a cooling instability for a hot gas since the gas 
would realistically have a broad range of the angular mo- 
mentum values. Additionally, a cloud with an initial size 
of a parsec or larger, tidally disrupted and shocked, should 
settle in a disk of a size comparable to its initial radius. 
There is of course a direct observational test which would 
distinguish between t he accretion dis k versus the compact 



infalling cloud idea oflSanderj l)1998D - one simply has to 
establish whether the orbits of the young stars in the two 
stellar rings are nearly circular or they are strongly eccen- 
tric. As we showed in 21 stars born in an accretion disk 
in Sgr A* would still retain their nearly circular orbits. 

The st a.ndard theory of gravitational i nstability for a 
thin disk l|Toomrel Il96l IPaczvhskil Il97i) predicts that 
the minimum mass of the gas in the disk that would 
make it gravitationally unstable for the parameters ap- 
propriate to our Galactic Center is ^ 10* Mq (0. It 
would be interesting to compare the predicted stellar 
mass resulting from star formation in such an accre- 
tion disk with the current stellar content of the rings. 



Unfortunately theoretical uncertainties for the efficiency 
of star formation in self-gravitating disks are too large. 
Ishi osman. Begelman. fc FrankI l|l99(t have shown that if 
the cooling time of a self-gravitating disk is shorter than 
tdyn — ^^^1 then the disk will fragment and form stars 
and/or planets. For longer coolin g times, it was argue d 
that the disk does not fragment ijShlosman et alJll99Cl|) . 
Numerical simu lations with a constant cooling time by 
confirmed this, and have shown that the 
disk settles into a stable state where the cooling is offset 
by the energ y input generate d by gravitational instabili- 
ties (see also |Paczvnskil[r978|) . Yet for disk temperatures 
of order 10^ Kel vin, the opacity is st rongly dependent on 
the temperature. I.Tohnson fc Gammie (2.0,0^) showed that 
in the non- linear stage of the instabilities, the local cooling 
time may be orders of magnitude smaller than that found 
in the unperturbed disk model. However AGN disks are 
usually hotter than this and hence the non-linear effect 
should be weaker. 

Nevertheless, we believe that Sgr A* accretion disk 
was likely consumed almost entirely in the star forma- 
tion episode rather than has been accreted by the SMBH. 
There is no doubt about star formation here: there are 
dozens of brightest and quite massive stars in each of the 
stellar rings with 3-D velocity measurements. There are 
additional numbers of dimmer stars that have only 2-D 
velocities measured but are strongly suspe cted of belong- 
ing to these same rings (IGenzel et al ] l2003i §3.7). The ac- 
cretion time scale on the embedded stars is very short (13 
and Figure |2l compared to the disk viscous time scale. 
We have also shown in this paper, that neither disk plane 
rotations, warps, or the stellar iV-body scattering (unless 
ther e is already more mass in th e stars than in the gas, see 
also' Cuadra fc Navakshinll2004^ can "shake" the stars off 
the disk midplane. In addition, each massive star opens up 
a radial gap in the accretion disk around it. These stars 
would not let the standard accretion disk to transfer the 
gas into the SMBH simply because they are in the way of 
the gas flow. 

While the standard accretion disk equations are not 
applicable to the region where the disk becomes self- 
gravitating, the stellar accretion time scale is shorter than 
the viscous time by 3-4 orders of magnitude. We experi- 
mented with a prescription for the accretion disk equations 
which introduces turbulent energy and pressure in addi- 
tion to the thermal ones to keep the disk marginally stable 
(i.e., Q ^ 1), and found that the turbulent energy content 
should be unrealistically high to reduce the accretion rate 
onto the stars sufficiently. 

We have also found (iQ that the geometrical thickness 
of the inner 2 — 4" stellar ring, and the stellar velocity dis- 
persion, could be explained by A'^body scatterings between 
the members of the same ring if the initial stellar mass was 
as high as 10* M(t,. However, the outer 4 — 7" projected 
distance ring l)Genzel et al.ll2003|) is too thick to result 
from the internal scatterings. We believe that a better ex- 
planation is stellar ring warping due to a non-spherical 
gravitational potential ([0, e.g. due to the presence of 
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the inner ring. Estimating the rate at which the rings get 
distorted, we tentatively set an upper limit on the time- 
average total mass of each of the gas-star disks (rings) at 
around 10^ AIq. Future numerical A'^body modeling and 
direct comparison to stellar orbits may tighten this limit. 

We hope that future observations of the stellar mass 
content in the two rings in Sgr A* , and also observations of 
the inner Galaxy ISM budget, could be used to constrain 
the initial mass of the gaseous accretion disk, and its fur- 
ther fate. As we have shown, the gaseous disk mass should 
have been in the range of (1 - 10) x IQ-^Mq, 10 to 100 
times higher than the present d ay mass in the observed 
stellar rings llCxenzel et alJl200,'^ . If most of the disk gas 
was used to make stars, then more of these stars and/or 
their remnants should be found in the future in the inner 
~ 0.2 pc of the Galaxy. In addition, one may look for ev- 
idence of a hot high metallicity bubble in the inner 1 kpc 
of the Galaxy produced by stellar winds and supernova 
explosions. 

If instead the gas was mostly accreted by Sgr A* , then 
there should be evidence of a past quasar phase. The re- 
quired accretion rate, (10'' - 10^) year = 10"^ 
to lO"'^ Mq year~^ is comparable with the Eddington ac- 
cretion rate for Sgr A*, MEdd ^ 0.03 Mq year^^. This 
would have to be a very rare event in Sgr A* recent life 
since the other nearby galaxies either have no AGN or 
have very weak ones with X-ray l uminosity usually sma ller 
than Lx ^ 10"^° erg sec"! Ce.g.. IZang fc Meur£ll200lh . A 
hot buoyant radio bubble would most likely be present 
in the Milky Way halo, as accretion onto the SMBH is 
widely believed to go hand in hand with superluminous 
jet outflows. 

10. Conclusions 

Our main results are as following: 

1. The minimum mass of each of the disks needed to 
form the observed young stars by self-gravity is around 
IO^Mq. 

2. The observed stellar velocity dispersions in the outer 
ring is too large to result from A'^body interactions be- 
tween stars belonging to the same ring. The orbital pre- 
cession of stars caused by the potential of the other disk 
can explain the observed disk thickness and velocity dis- 
persion if the time average stellar and gaseous mass in the 
inner disk is in the range (3 — 10) x 10'* Mq. 

3. Few million years ago, Sgr A* had a good chance to 
become a very bright AGN with the bolometric luminosity 
L ~ 10^"* — 10"*^ erg/sec, but was robbed of most of its 
gaseous fuel by nuclear star formation in a self-gravitating 
accretion disk. Nevertheless, even if only a few percent of 
the available disk fuel was captured by Sgr A* , the SMBH 
in our Galactic Center was as bright as i - lO''^ - 10^^ 
erg/sec. 

We have also shown that capture of stars from the 
"old" relaxed isotr opic Sgr A* star cluster (the cusp; see 
iGenzel et al.ll200^ is inefficient unless the gaseous disk 



mass were as high as 10^ Mq. In addition, the role of pos- 
sible accretion disk midplane changes was estimated. It 
was found that the embedded stars inertia would have 
been efficient in taking the stars out of the body of the 
disks only if the disk plane changes its orientation on time 
scales much shorter than the disk viscous time. 

Acknowledgements. The authors acknowledge useful discus- 
sions with Eugene Churazov, Walter Dehnen, Reinhard Genzel, 
Marat Gilfanov, Peter Predehl and Rashid Sunyaev. 



References 

Artymowicz, P. 1994, ApJ, 423, 581 

Artymowicz, P., Lin, D. N. C., & Wampler, E. J. 1993, ApJ, 
409, 592 

Bate, M. R., Lubow, S. H., Ogilvie, G. I., & Miller, K. A. 2003, 

MNRAS, 341, 213 
Binney, J. & Tremaine, S. 1987, Galactic dynamics (Princeton 

University Press) 
Chandrasekhar, S. 1943, ApJ, 97, 255 
Collin, S. & Zahn, J. 1999, A&A, 344, 433 
Cuadra, J. & Nayakshin, S. 2004, in Growing Black Holes: 

Accretion in a Cosmological Context, ed. A. Merloni, S. 

Nayakshin & R. Sunyaev (in press) 
Fryer, C. L., Woosley, S. E., & Heger, A. 2001, ApJ, 550, 372 
Gammie, C. F. 2001, ApJ, 553, 174 

Genzel, R., Schodel, R., Ott, T., et al. 2003, ApJ, 594, 812 
Gerhard, O. 2001, ApJ, 546, L39 

Chez, A. M., Duchene, G., Matthews, K., et al. 2003, ApJ, 

586, L127 
Goodman, J. 2003, MNRAS, 339, 937 
Goodman, J. & Tan, J. C. 2004, ApJ, 608, 108 
Hansen, B. M. S. & Milosavljevic, M. 2003, ApJ, 593, L77 
Johnson, B. M. & Gammie, C. F. 2003, ApJ, 597, 131 
Kim, S. S., Figer, D. F., & Morris, M. 2004, ApJ, 607, L123 
Kim, S. S. & Morris, M. 2003, ApJ, 597, 312 
Kolykhalov, P. I. & Sunyaev, R. A. 1980, Soviet Astron. Lett., 

6, 357 

Koyama, K., Maeda, Y., Sonobe, T., et al. 1996, PAS J, 48, 249 
Krabbe, A., Genzel, R., Eckart, A., et al. 1995, ApJ, 447, L95 
Kulsrud, R. M., Mark, J. W. K., & Caruso, A. 1971, Ap&SS, 
14, 52 

Levin, Y. & Beloborodov, A. M. 2003, ApJ, 590, L33 
Lissauer, J. J. 1987, Icarus, 69, 249 
Liszt, H. S. 2003, A&A, 408, 1009 

Maillard, J. P., Paumard, T., Stolovy, S. R., & Rigaut, F. 2004, 
A&A, 423, 155 

McMillan, S. L. W. & Portegies Zwart, S. F. 2003, ApJ, 596, 
314 

Milosavljevic, M. & Loeb, A. 2004, ApJ, 604, L45 

Morris, M., Ghez, A. M., & Becklin, E. E. 1999, Advances in 

Space Research, 23, 959 
Nayakshin, S. 2004, MNRAS, 352, 1028 

Nayakshin, S., Cuadra, J., & Sunyaev, R. 2004, A&A, 413, 173 
Ostriker, J. P. 1983, ApJ, 273, 99 
Paczyhski, B. 1978, Acta Astron., 28, 91 

Rees, M. J. 2002, in Lighthouses of the universe: the most 
luminous celestial objects and their use for cosmology, ed. 
M. Gilfanov, R. Sunyaev & E. Churazov (Berlin: Springer), 
345 



10 



S. Nayakshin and J. Cuadra: Sgr A* few Million years ago 



Revnivtsev, M. G., Churazov, E. M., Sazonov, S. Y., et al. 
2004, ArXiv Astrophysics e-prints, astro-ph/0408190 

Sanders, R. H. 1998, MNRAS, 294, 35 

Schodel, R., Ott, T., Genzel, R., et al. 2002, Nature, 419, 694 
Shakura, N. I. & Sunyaev, R. A. 1973, A&A, 24, 337 
Shlosman, I. & Begelman, M. C. 1989, ApJ, 341, 685 
Shlosman, I., Bogolman, M. C, & Frank, J. 1990, Nature, 345, 
679 

Sunyaev, R. A., Markevitch, M., & Pavlinsky, M. 1993, ApJ, 
407, 606 

Svensson, R. & Zdziarski, A. A. 1994, ApJ, 436, 599 
Syer, D., Clarke, C. J., & Rees, M. J. 1991, MNRAS, 250, 505 
Tanga, P., Weidenschilling, S. J., Michel, P., & Richardson, 
D. C. 2004, ArXiv Astrophysics e-prints, astro-ph/0408441 
Toomre, A. 1964, ApJ, 139, 1217 
Yu, Q. & Tremaine, S. 2002, MNRAS, 335, 965 
Zang, Z. & Meurs, E. J. A. 2001, ApJ, 556, 24 



